Radiative and Collisional Jet Energy Loss in a Quark-Gluon Plasma 



G.Y. Qin 1 ^ J- Ruppert 1 ' 3 , C. Gale 1 , S. Jeon 1 , G. D. Moore 1 , and M. G. Mustafa 2 

1 Department of Physics, McGill University, Montreal, QC, USA 2T8, Canada 

2 Saha Institute of Nuclear Physics, 1/AF Bidhannagar, Kolkata, 700064, India 

3 Institut fur Theoretische Physik, J. W. Goethe University Frankfurt, 
Max-von-Laue-Str. 1, D-60438 Frankfurt am Main, Germany 



We calculate radiative and collisional energy loss of hard partons traversing the quark- 
gluon plasma created at RHIC and compare the respective size of these contributions. We 
employ the AMY formalism for radiative energy loss and include additionally energy loss 
Q\ ' by elastic collisions. Our treatment of both processes is complete at leading order in the 

fN| , coupling, and accounts for the probabilistic nature of jet energy loss. We find that a solution 

of the Fokker-Planck equation for the probability density distributions of partons is necessary 
for a complete calculation of the nuclear modification factor Raa for pion production in heavy 
ion collisions. It is found that the magnitude of Raa is sensitive to the inclusion of both 
collisional and radiative energy loss, while the average energy is less affected by the addition 
|' of collisional contributions. We present a calculation of Raa for ir° at RHIC, combining our 

energy loss formalism with a relativistic (3+l)-dimensional hydrodynamic description of the 
. thermalized medium. 



I. INTRODUCTION 
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■<^j- ■ The main goal of relativistic heavy ion collisions is to study strongly interacting nuclear matter 

at extreme temperatures and densities and investigate the possibility of a phase transition between 
hadronic matter and a quark-gluon plasma (QGP). Experiments at the Relativistic Heavy Ion 
Collider (RHIC) have demonstrated a significant suppression of high pr hadrons produced in 
central A + A collisions in comparison to those in binary scaled p + p collisions [j], 0]. This 
result has been referred to as jet quenching and has been attributed to the energy loss of hard px 
partons in a hot quark-gluon plasma [3[. Experimentally, jet quenching can be characterized by 
' measurements of various quantities, i.e., the nuclear modification factor Raa, the elliptic flow V2 

at high pt, and high px hadron correlations, etc. 

In [J] we presented a systematic calculation of Raa for 7r° in central and non-central Au+Au 
collisions at y/s — 200 AGeV combining the Arnold, Moore and Yaffe (AMY) [f| formalism with 
a (3+l)-dimensional ideal hydrodynamical model [6]. In this contribution, we consistently incor- 
porate both collisional and radiative energy loss in the same framework [tJ ■ 



II. JET ENERGY LOSS IN THE AMY FORMALISM 



In the AMY formalism, the evolution of the jet momentum distribution P(E,t) — dN(E,t)/dE 
in the QGP is obtained by solving a set of coupled Fokker-Planck type rate equations for quarks 
plus anti-quarks and gluons, which have the following generic form, 
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where dT(E,ui) / 1 dudi is the transition rate for processes in which partons of energy E lose energy 
u). The <jj < part of the integration incorporates processes which increase a particle's energy. 
The radiative part of the transition rates is discussed in 0, HI, [9[ . In the following, we address the 
calculation of elastic collision rates Q. 

As collisional energy losses are more sensitive to small energy transfers in comparison to radiative 
ones, it should be an adequate procedure to approximate elastic collisions by a mean energy loss 
plus a momentum diffusion term as dictated by detailed balance. If elastic collisions turn out to 
dominate jet quenching we may want to improve this treatment, but if they are subdominant it 
should be sufficient to quantify their effect. 

As we perform the numerical computation of discretized Eq. |T]) such that J duj — > A w X)w=«A > 
the transition rates for elastic collisions are approximated by two spikes at ui = ±A W , 

with Jb the Bose-Einstein thermal distribution function. For small A w , the above collisional 
transition rates yield the correct energy loss rate and preserve the detailed balance. In the small 
5u limit, this procedure to incorporate collisional energy loss is equivalent to introducing the drag 
term (dE/dt)dP(E)/dE, and the diffusion term T(dE/dt)d 2 P(E)/dE 2 into Eq. (TJJ. 

At leading order, the mean energy loss by elastic collisions is calculated in kinetic theory, 
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where C co ii and D co \\ are 0(1) constants depending on different channels Q, and m 2 = 2na s T 2 (l- 
Nf/6) is the thermal gluon mass. 



III. RESULTS AT RHIC 

Now we employ this formalism to a realistic description of energy loss of hard pt leading partons 
in the soft nuclear medium in Au+Au collisions at y/s = 200 AGeV at RHIC. We utilize a fully 
(3+l)-dimcnsional hydrodynamical model for the description of the thermalized medium created 
in the collisions as it has been shown to give a good description of bulk properties at RHIC @ . 
The product of the initial hard parton densities is determined from the overlap geometry between 
two nuclei in the transverse plane of the collision zone. The initial momentum distribution of hard 
jets is computed from perturbative QCD, using the factorization theorem. The evolution of jet 
momentum distribution in the thermalized medium is evaluated in AMY [5|, solving Eq. fTJ with 
both collisional and radiative energy loss. The final hadron spectrum at high pt is obtained by the 
fragmentation of jets in the vacuum after their passing through the (3+l)-dimensional expanding 
medium. For further details see Q where the radiative energy loss has been studied in an analogous 
manner. 

The nuclear modification factor Raa is defined as the ratio of the hadron yield in A+A collisions 
to that in p+p interactions scaled by the number of binary collisions N co vC- 

r,h / j — ! 1 dN h (b)/d 2 p T dy 

R h AA (b,p T ,y) = AT , M , AT l " —, ■ 4 

Acoii(o) dN£ p /d z p T dy 

In the AMY formalism, the strong coupling constant a s is the only quantity which is not uniquely 
determined by the model, once the temperature and flow profiles are fixed by the initial conditions 
and subsequent evolution of the (3+l)-dimcnsional hydrodynamics. In this work, we take its value 
to be constant at a s — 0.27, which reproduces the most central data (see Fig. [2]). 




FIG. 1: (Color online) A single quark jet starts from the center and propagates in plane through the 
nuclear medium created in most central Au+Au collisions (b — 2.4 fm) at RHIC. Left: The final probability 
distribution P(E,tf) of a single quark with initial energy Ei — 16 GeV after passing through the medium. 
Right (from [7(): The fractional averaged energy loss of a quark jet with initial energy Ei after passing 
through the medium. 



To illustrate how collisional and radiative energy loss influence the time evolution of the leading 
parton distributions, we first consider an initial single quark jet starting from center and propagat- 
ing in plane through the thermalized medium created in central collisions (b — 2.4 fm) at RHIC, 
as shown in Fig. [1] In the left panel, we compare final parton distribution P(E,tf) under three 
different situations: (1) with only collisional energy loss, (2) with only radiative energy loss, and 
(3) with both energy loss mechanisms. In the right panel, the mean energy loss of quark jets after 
passing through the medium is shown as a function of their initial energy Ei . As expected, induced 
gluon radiation leads to a much larger mean energy loss than elastic collisions, in agreement with 
fig ]. While the average energy is not very affected by the addition of collisional contributions, large 
differences are observed for the time evolutions of P(E,t) between case (3) and case (2). This is 
especially true for energies close to the initial parton energy Ei . As the initial parton spectrum in 
relativistic nucleus-nucleus collisions is steeply falling, stronger differences in Raa can result. 

In the left panel of Fig. [21 we present the calculation of Raa for 7r° measured at mid-rapidity 
for two different impact parameters, 2.4 fm and 7.5 fm, compared with PHENIX data for the 
most central (0-5%) and mid-central (20-30%) collisions Q. We present Raa for purely collisional 
(1) and purely radiative (2) energy loss, as well as the combined case (3). One finds that while 
the shape of Raa for case (3) is not strongly different from case (2), the overall magnitude of 
Raa changes significantly. The magnitude of Raa is therefore sensitive to the actual form of the 
parton distribution functions at fragmentation and not only to the average energy loss of single 
partons (compare Fig. [I]). In 4], a systematic study of the observational implications on Raa 
measurements due to only radiative energy loss was presented. Here we recalculate Raa versus 
reaction plane including elastic collisions, as shown in the right panel of Fig. [5J and only small 
differences (after rescaling the coupling strength from a s — 0.33 to a s = 0.27) is found in the 
shape of Raa for the pr range explored. 



IV. SUMMARY 



In conclusion, we have calculated radiative and collisional energy loss of hard partons in the hot 
and dense nuclear medium being created at RHIC. The induced gluon bremsstrahlung is treated 
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FIG. 2: (Color online) Left (from 0|, 0]) : Raa for 7r° in most central and mid-central Au+Au collisions 
at RHIC. The dashed curves includes radiative energy loss only, the dash-dotted curves include collisional 
energy loss and the solid curves include both radiative and collisional energy loss. The dotted curves are 
the results from |4[ with radiative energy loss only (a a = 0.33). Right: Raa for tt in plane and out of 
plane in Au+Au collisions (b = 7.5 fm) at RHIC before and after the inclusion of collisional energy loss. 



in the AMY formalism and elastic collisions are incorporated using a drag plus diffusion term 
reproducing average energy loss rate and detailed balance. While the additional average energy 
loss induced by elastic collisions is small in comparison to the radiative one, the time evolutions of 
the parton distributions in both cases differ significantly. We find that the inclusion of collisional 
energy loss significantly influences the overall magnitude of Raa for tt° at RHIC, while the shape 
of Raa does not show a strong sensitivity. 
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